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We examine profiles and scaling properties of the entropy of the intergalactic gas in 
a sample of 66 virialized systems, ranging in mass from single elliptical galaxies to rich 
clusters, for which we have resolved X-ray temperature profiles. Some of the properties 
we derive appear to be inconsistent with any of the models put forward to explain the 
breaking of self-similarity in the baryon content of clusters. In particular, the entropy 
profiles, scaled to the virial radius, are broadly similar in form across the sample, 
apart from a normalization factor which differs from the simple self-similar scaling 
with temperature. Low mass systems do not show the large isentropic cores predicted 
by preheating models , and t he high entropy excesses reported at large radii in groups 
bv lFinoeuenov et al.1 l)2002|) are confirmed, and found to extend even to moderately 
rich clusters. We discuss the implications of these results for the evolutionary history 
of the hot gas in clusters, and suggest that preheating may affect the entropy of 
intracluster gas primarily by reducing the density of material accreting into groups 
and clusters along cosmic filaments. 
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1 INTRODUCTION 

In the widely accepted standard model for cosmic structure 
formation, the Universe evolves hierarchically, as primor- 
dial density fluctuations, amplified by gravity, collapse and 
merge to form progressively larger systems. This hierarchical 
development leads to the prediction of self-similar scalings 
between systems of different masses and at different epochs. 
These scalings are also seen in cosmological simulations in- 
volving only gravitationally driven evolution, including com- 
pression and s hock heating of the baryonic matter. Su ch sim- 
ulations (e.g. iNavarro et allll995l: iFrenk et al]ll999l) result 
in haloes in which the density profiles of both dark mat- 
ter and baryonic material, when radially scaled to the virial 
radius (which we define here to be the radius ii200, within 
which the mean density of a system is 200 times the criti- 
cal density of the Universe) are almost identical in virialized 
systems covering a wide range of masses, from individual 
galaxies to rich clusters. 

Given self-similar scalings of gas temperature and den- 
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sity, scaled X-ray surface brightness profiles are also ex- 
pected to be similar. Furthermore, a simple scaling is 
expected between X-ray luminosity, Lx, and tempera- 
ture T. Assuming that the emission is dominated by 
bremsstrahlung, Lx <x M^^R'^^T^^^, or Lx oc fgasT^t 
where Mgas is the gas mass within radius R, and fgas = 
Mgas/M is the gas mass fraction. X-ray properties of clus- 
ters deviate substantially f rom this s imple scalin g, and 
the o bserved Lx -T relation jWhite et aL.1997; .Markevitcbl 
Il998ll is considerably s teeper than in the clus ter regime, 
and ste epens further j Hglgd on fc Ponmanll200nh in galaxy 
groups. iPonman et al showed that the latter effect 

is due to the suppression of Lx in galaxy groups, arising 
from a reduction in gas density in the inner regions of poor 
systems, relative to richer clusters. 

It is instructive to view this in terms of the entropy of 
the intergalactic medium (IGM), which in the self-similar 
case should increase in a very simple scaling with the mean 
temperature of viria l ized systems. In practice , it is found 
l|Ponman et alJll999l: iLlovd-Davies et aljEoonl) that an ex- 
cess in the entropy, above the self-similar prediction, is ap- 
parent in the inner regions of poor clusters and groups, out- 
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side the dense central regions in which cooling is expected to 
radiate away entropy within the age of the Universe. This 
effect has been referred to as the 'entropy floor', with the 
implication that additional physical processes, beyond grav- 
ity and resulting compression and shock heating, have acted 
to set a lower limit to the entropy which the gas in collapsed 
haloes can have. 

A great deal of theoretical work has been devoted to ex- 
plaining this phenomenon over the past few years. As we will 
discuss in some detail later, the explanations proposed fall 
into three main classes: the gas heis been heate d either at an 
early epoch, before clusters were assembled |^^igei[^^l 



Evrard & Henrv "l991; 'Cavalicre et al.' 'l997; 'Balo gh et all 
'l99a:.Valageas fc Silk 1999: Tozzi fc Norman.200L) . or it has 



been heated in situ by star forma tion and /or energy input 
from active galactic nuclei (AGN) ('Bowei'1997'; Loewenstein 
I2OOO.; Vgit & Brvan 2001; Nath & Rovchowdhurv 2002). Al- 
ternatively, some authors have argu ed (|Kni£liO^.^Qnmani 
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Wu fc Xuell2002al : iDave et al^l2002^ that cooling alone will 
remove low entropy gas from the centres of haloes, producing 
a very similar effect to non-gravitional heating. 

In the present paper we aim to confront these models 
with the observed properties of the hot gas in a large sample 
of galaxy systems, spanning a wide range in total mass. In 
the fullness of time, high quality X-ray observations of the 
density and temperature structure of the IGM will be avail- 
able from XMM-Newton and Chandra. However, at present, 
such observations are sparse, and it is essential to have a 
broadly representative and wide-ranging sample of virial- 
ized systems in order to study scaling properties. The value 
of this in the context of similarity-breaking in clusters has 
already been shown by a number of earlier studies. 

In the present, and companion papers. 'Sanders on et alJ 
(1200,1) (Paper I) and Sanderson & Ponman ( 2003) (Pa- 
per II), we examine scaling properties derived from the 
largest sample of virialized systems with resolved X-ray tem- 
perature profiles yet assembled. Following a brief description 
of the sample and our analysis in section |5| (details are given 
in Paper I), we present the profiles and scaling properties 
for the entropy and temperature across our sample in sec- 
tions|3and|l| These results are used in section|21 along with 
relevant results from Papers I and II, to test the various 
models proposed to account for the entropy fioor, and fi- 
nally in section |S| we draw our conclusions from this study, 
and propose a new model to explain the behaviour of the 
entropy of intracluster gas. 



2 SAMPLE AND ANALYSIS 

Our sample comprises 66 virialized systems, from rich clus- 
ters of galaxies, through groups a nd down to the l e vel of 
individua l galaxy- sized haloes. In ISanderson et alj i2003l 
hereafter iPaoer l|l . we reported a detailed study of the 
3-dimensional X-ray properties of this sample, based on 
data from the ROSAT and ASCA observatories, which we 
assembled from the work of th ree separate investigators 
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combined with a number of cool groups analysed specially 
to provide better coverage of the crucial low end of the mass 
range. To each system we fitted analytical functions, describ- 
ing the gas density and temperature variation with radius, 
outside any central cooling region, which was excised or fit- 
ted separately. This approach allows us to put the X-ray 
data from the three earlier studies on a unified footing, and 
gives us the freedom to extrapolate the gas properties to 
arbitrary radius. We used a beta model to parametrize the 
density and specified the temperature variation with either a 
linear ramp or a polytropic IGM. We have used these data to 
determine the gravitating mass profile and thus to calculate 
radii of overdensity in a self-consistent manner. Similarly, 
we have derived mean temperatures for each system, by av- 
eraging the gas temp erature within 0.3-R200, weighted by its 
luminosity 



: gas temp i 
jPaper J) . 



3 ENTROPY AND TEMPERATURE 
DISTRIBUTIONS 

For convenience, we define 'entropy' in the present paper as 



S = T/nl^^ keV cm 



(1) 



which relates directly to observations. This has been referred 
to by a number of authors as the 'adiabat', since (apart from 
a constant relating to mean particle mass) it is the coefficient 
relating pressure and density in the adiabatic relationship 
P = Kp''. Hence S is conserved in any adiabatic process. 
Note that the true thermodynamic entropy is related to our 
definition via a logarithm and additive constant. 

In Fig.0 we overlay the scaled entropy profiles for all 66 
systems in our sample. Under the assumption that all these 
systems form at the same redshift, their mean mass densities 
should be identical. Hence in the simple self-similar case, 
where all have similar profiles and identical gas fractions, S 
will simply scale with mean system temperature T. We apply 
this scaling, and scale the radial coordinate to -R200 for each 
system, derived from our fitted models (Paper I). It can be 
seen that the entropy profiles of the cooler systems, scaled in 
this way, tend to be significantly higher than those of richer 
clusters. To see this more clearly, in the right panel of Fig.0 
we show the profiles grouped into bands of mean system 
temperature. In this grouped plot we have excluded the two 
galaxies in our sample, whose profiles may be dominated by 
stellar wind losses, rather than by the processes operating 
in groups and clusters. 

It can be seen that there is a strong tendency for the 
scaled entropy to be higher, at a given scaled radius, in cooler 
systems. Simulations and analytical models of cluster forma- 
tion involving only gravity and shock heating, produce en- 
t ropy pro files with logarithmic slopes of approximately 1.1 
frozzi XNormanll20"0ll) . which agrees rather well with the 
slope of our profiles outside 0.17?2oo- The mean profile in the 
2.9-4.6 keV band is significantly affected by two systems, 
visible in the left hand panel, which have strongly rising 
profiles. One of these is AWM7, for which the temperature 
profile is subject to especially large systematic uncertain- 
ties, as a result of the strong cooling fiow, as discussed in 
P aper I. The temperature m odel adopted, from the analysis 
of iLlovd-Davies et al.l i2000l) . rises strongly with radius. As 



Entropy and similarity in galaxy systems 3 




Figure 1. Entropy profiles for each system in our sample (left panel) derived from our fitted models, each scaled by 1/T. The line style 
of each profile denotes the mean temperature of the system, as described below. In the right panel, these profiles (excluding the two 
galaxies, as discussed in the text) have been grouped into temperature bands: 0.3-1.3 keV (solid), 1.3-2.9 keV (dashed), 2.9—4.6 keV 
(dotted), 4.6-8 keV (dot-dash) and 8-17 keV (dot-dot-dot-dash). The bottom line shows the slope of 1.1 expected from shock heating. 
Its normalization is arbitrary. 



a result, the slope of S{r) in this band is almost certainly 
overestimated in the figure. 

One possibility which we can discard right away, is that 
the observed scaling results from systematic differences in 
the formation epochs of low and high mass systems. In hier- 
archical models, low mass systems are expected to form, on 
average, earlier than high mass ones. Hence they should tend 
to have higher mean densities, and therefore lower gas en- 
tropy - the opposite of what we observe. It is true that, due 
to the usual selection effects, the rarer and more luminous 
massive systems in our sample tend to be at higher redshifts 
than the groups. Hence if, contrary to expectations, all sys- 
tems virialized at the redshift we observe them, then the 
more massive systems would be more dense, and have lower 
entropies. However, this effect is considerably smaller than 
the trend which we observe. The highest redshift clusters in 
our sample are at z ~ 0.2, hence their mean densities, scal- 
ing as [l + Zf)^, could be 70% higher than the lowest redshift 
systems, and their entropies (scaling as n^^^) consequently 
lower by a factor of 1.44. All but four members of our sam- 
ple lie at z < 0.1, and hence would have entropies changed 
by less than 20% as a result of such a redshift- dependent 
density scaling. In contrast, Fig.^shows that the scaled en- 
tropy actually differs by a factor 3 between the high and 
low temperature bins for our sample. 

At small radii, our fitted models exclude the effects of 
cooled gas, since any central cooling region is excised, or rep- 
resented by a separate component in our analysis (secLj^J. 
This central cooling region is present in 54 of our 66 sys- 
tems (see Table 1 in Paper I), and in these cases its ra- 
dius (rcooi) ranges from 0.03-R200 to 0.2_R2oo, with a me- 
dian value of Tcooi = 0.06-R200- Apart from the coolest sys- 
tems, our entropy profiles generally fiatten inside 0.1-R200- 
This effect i s also seen in many cosmological simulations 
jFrenk et al. 1999) even in the absence of non-gravitational 
heating and cooling processes, and appears to result from 
the introduction of a core into the gas density distribution, 
due to transfer of energy between baryonic and dark matter 



during merger events jEke et alJll998l) . Preheating models 
generally predict large isentropic central regions in low mass 
systems, which are not seen in our data. We will return to 

this point in section |S| below. 

The preheating model of iDos Santos fc Dorl (|20o3) 
predicts that entropy should scale according to (1 -I- T/Tq), 
rather than T, where To is a constant related to the degree 
of preheating, which they estimate as T q=2 keV to provide 
a best fit to the entropy floor data of iLlovd -Davies et alJ 
i2000l) . In Fig. |^ we show the effect of this scaling on our 
temperature-grouped entropy proflles. This scaling does in- 
deed bring the profiles into good agreement, apart from the 
fact that our coolest systems show little sign of any central 
entropy core. Note from the left hand panel in Fig. Q that 
this is a general feature of almost all the entropy profiles for 
cool systems, rather than a result of averaging together sys- 
tems with cores and others with strongly dropping central 
entropies. 

It is also instructive to compare temperature profiles 
across the range of masses in our sample. We use the virial 
radii and masses calculated from our fitted models to com- 
pute the virial temperature 



r2oo = 



GMnTUp 

2ii200 



(2) 



expressed in keV. This is used to normalize each tempera- 
ture profile, and the scaled profiles are grouped in tempera- 
ture bands, to reduce the large scatter and make trends more 
obvious. The result is shown in Fig-O] As with the grouped 
entropy plots, we have excluded the two galaxies from these 
proflles. We have also omitted AWM7 from the 2.9-4.6 keV 
band, since its highly uncertain (see discussion above) and 
strongly rising T(r) proflle distorts the mean profile for the 
whole band. 

Raising the entropy of the IGM results in increased tem- 
peratures as well as lower densities, and for a given level of 
entropy increase this effect will be most prominent in low 
mass systems, where the 'natural' shock-generated entropy 
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Scaled Radius (R/R200) 



Figure 2. The variation of gas entropy (scaled by (1 + T/Tq)"^) witli scaled radius, grouped by system temperature. The solid line 
represents the coolest systems (excluding the two galaxies) (0.3-1.3 keV), increasing in temperature through dashed (1.3—2.9 keV), dotted 
(2.9-4.6 keV), dot-dashed (4.6—8 keV) and finally dot-dot-dot-dashed (8—17 keV). The lower solid line (with arbitrary normalization) 
indicates the slope of 1.1 expected from shock heating. 



is lower. As can be seen, our observations do indeed show 
such an effect at large radii - at i?200 the scaled tempera- 
tures are ordered in precisely this way. However, at smaller 
radii, this is not the case, and in particular, the coolest sys- 
tems display very flat temperature profiles, whilst the most 
massive clusters show the largest rise in temperature be- 
tween i?200 and the centre. This behaviour is not what is 
predicted by most models, as we will discuss later. 



4 SCALING PROPERTIES 

The c laim of the discovery of an entropy floor in galaxy sys- 
tems jPonman et alJ^1999^ was based on the measurement 
of gas entropy at a scaled radius of O.I-R200, in systems span- 
ning a wide temperature range. This radius was chosen to 
lie close to the centre, where shock-generated entropy should 
be a minimum, hence maximising the sensitivity to any ad- 
ditional entropy, whilst lying outside the region where the 
coohng time is less than the age of the Universe, and hence 
th e entropy may be reduced. T his initial study was improved 
by iLlovd-Davies et alj l|20 0(]|). who avoided the isothermal 
assumption made bv lPoiinmnet''all and derived an entropy 
floor value of 139 h^o^^ keV cm^ from a sample of 20 sys- 
tems, which is essentially a subset of the present work. 

In Fig. 21 we show the corresponding plot from our much 
larger sample. With the beneflt of this increase in sample 
si ze, the trend looks rather different from its appearance 
in lLlovd-Davies et alJ 12000.,') . where a dearth of systems in 
the 1.5-3.5 keV band led to the interpretation of a relation 
which followed the self-similar line down from hot systems, 
flattening towards a floor value of 139 h^Q^^ keV cm^, cor- 
responding, with our value of Ho, to 124/170"^/^ keV cm . 



However, it appears from our data that the behaviour is 
rather that of a slope in S{T) which is significantly shal- 
lower than the self-similar relation S oc T throughout. Us- 
ing unweighted orthogonal regression lllsobe et aljil990l) . we 
obtain a relat ion g(0. 1^200) oc r O-65±o.05^ which is marked 
in Fig.|l| The lLlovd-Davies et alJ floor value lies at the bot- 
tom of this trend, but it is not clear whether it sets a lower 
limit to the entropy in galaxy groups, since no groups with 
T<0.6 keV have been bright enough for detailed study to 
date. The two galaxies in our sample do clearly have en- 
tropies which lie below the floor level, however much or 
all of the hot gas in these systems may have its origin in 
stellar m ass loss, rathe r than retained primordial material 
f O ' Sulliv an et ai] |2001^ , so it may well be fortuitous that 
they fall close to the trend set by the groups and clusters. 

Grouping the points into temperature bins (Fig. |SJ 
makes the unbroken nature of the trend very clear. An un- 
weighted orthogonal flt to these grouped points gives a loga- 
rithmic slope of 0.57±0.04, slightly flatter than that derived 
from Fig. |1| Note that these results suggest the presence of 
excess entropy of ~ 100 keV cm^ in all temperature bands, 
relative to the hottest systems. Fig. |S] also compares the 
effect on the relation of calculating -R200 from a measure- 
ment of the mass profile, compared to assuming a scaling 
R200 oc T^/^, as w as done (employing t he f ormula derived 
from s imul ations bvlNava rro et all ^1995^ 'l bv lPonman et all 
(I99J) a nd iLlovd-Davies et alj ll200(i) ."As discussed in Pa- 
per I, the lNavar^^^ilJ formula agrees reasonably well with 
our measurements in rich clusters, but in cool systems, the 
upward bias in temperature, r elative to simulations such as 
those of iNavarro et al.1 (|^93), which include only gravita- 
tional physics, causes the 7?200 oc T^/'^ scaling to overesti- 
mate R200 leading in turn to an upward bias in the entropy 
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Figure 3. The variation of gas temperature (scaled by T200) with scaled radius, grouped by system temperature. The solid line 
represents the coolest systems (excluding the two galaxies) (0.3—1.3 keV), increasing in temperature through dashed (1.3—2.9 keV), 
dotted (2.9-4.6 keV, excluding AWM 7), dot-dashed (4.6-8 keV) and finally dot-dot-dot-dashed (8-17 keV). 



(since 5* rises with radius). This may have contributed, in 
previous studies, to the appearance of flattening in the ^(r) 
relation towards low T. 

Whilst measurements close to the cluster centre provide 
the most sensitive probe of excess entropy, detection of ad- 
ditional entropy at large radii is especially interesting, since 
many preheating models predict that the rise in entropy is 
essentially restricted to those central regions where shock- 
generated ent ropy is le ss than the floor value set by preheat- 
ing. Finogu enov et alJ y002) were the first to find evidence 
for excess entropy at a much larger radius, R500, correspond- 
ing (Paper II) to ~ |7i200- They argued that the high excess 
entropy observed at large radii in groups and poor clusters 
indicates that their IGM is dominated throughout by the 
effects of preheating, with shocks playing little or no role. 
We aim, in this paper, to check this result with our larger 
sample. 

Our data, shown in Fig.|^ confirm the existence of sub- 
stantial excess entropy at i?5oo, above a self-similar extrap- 
olation from the values seen in rich clusters (dotted line). 
The trend is more clearly seen in temperature-grouped data 
shown in the right hand panel, and the logarithmic slope 
of S{T) is similar to that seen at 0.1i?200- The departure 
from self-similar entropy values is apparent not only in the 
coolest systems, but also in quite rich clusters. In fact, the 
largest absolute values of excess entropy (~ 1000 keV cm^) 
are seen in clusters with T ~3-4 keV. 

For direct comparison with iFinoeuenov et all i2002l) . 
we show in Fig.|7|the entropy at i?5oo scaled by M^^^^^ , and 
grouped into mass bins to suppress fiuctuations. For a set 
of self-similar systems, virialising at the same epoch (and 
hence having the same mean density), this scaling should 
renormalize the entropy to a constant value, independent 



of system temperature (since S (xT oc M^/^). Clearly real 
clusters do not follow this self-similar law. Whilst our plot is 
broad ly consistent in shape with Fig. 3 of lFinoeuenov et alJ 
i2002h . o ur larger sample again re veals important additional 
features. IFinoeuenov et alJ (l2002l) concluded that excess en- 
tropy was only present in systems with M500 < 10^"* K^^Mq 
(7 X 10^^ Mq for our choice of Hq), whereas it is clear from 
our data that a trend in scaled entropy is present across the 
full mass range of clusters and groups. 

IFinoeuenov et al] ll2002ll also examined the entropy at 
a fixed value of enclosed gravitating mass (3x 10" H-^Mq) 
over their sample of systems. The infall velocity of gas into 
an accretion shock should be similar to the free fall velocity 
at the shock radius, which depends upon the enclosed mass 
and mean density of the Universe at the epoch in question. 
In clusters, an enclosed mass of 3 X 10^^ h-^oMQ lies deep 
within the system, whilst in small groups it can lie close to 
i?200- Since cluster cores are generally assembled at higher 
redshifts than the outskirts of groups, one expects the gas 
density in the shell under consideration to be higher, and the 
entropy generated by the thermalised kinetic energy is there- 
fore lower, by a factor {1 -\- Zf), where is the redshift at 
which the shell was accreted. In practice, IFinoeuenov et aD 
I2OO23 found the entropy of this shell to be bimodal across 
their sample, with a typical value of ~300 keV cm^ for sys- 
tems with T ^ 3 keV, and with considerably lower values 
(scattered over the range 100-300 keV cm^) in hotter sys- 
tems. Such a distribution c annot be accounte d for on the 
basis of a 1-1-0/ scaling, and IFinoeuenov et al.l sueeested in- 
stead that the entropy in cool systems is set by preheating, 
taking place at z ^ 3, after many cluster cores had already 
collapsed. 

The corresponding plot for our sample, which is almost 



6 T. J. Ponman, A. J. R. Sanderson and A. Finoguenov 



o 
o 
o 



> 

CD 

>^ 
Q. 
O 



0) 
D 
O 




Temperature (keV) 

Figure 4. Gas entropy at O.I-R200 a function of system temperature. Tlie diamonds represent tlie two galaxies. The solid line is the 
best fit to the barred points. The dotted line has a self-simil ar slope of 1 and is normalized to the mean of the hottest 8 clusters. The 
dashed line is the entropy floor of 124 /i7o~^''^ keV cm^ from lLlovd-Davies et al 
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Figure 5. Gas entropy at O.I-R200 as a function of system temperature, excluding the two galaxies and grouped to a minimum of 8 
points per bin. Barred crosses are bas ed our measured valu es of -R20O1 whilst diamonds result from calculation of -R200 from mean system 
temperature using the T^^'^ scaling of lNavarro et all ^^9^. The dotted line shows the self-similar slope of 1, normalized to the 8 hottest 
clusters. 
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Figure 6. Entropy at -R500 as a function of system temperature for each of tlie 66 systems, with the galaxies marked as diamonds (left 
panel) and grouped to a minimum of 8 points per bin, excluding the galaxies (right panel). The dotted line shows the self-similar slope 
of 1, normalized to the 8 hottest clusters. 
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Figure 7. Gas entropy at -R50O1 normalized by Afgpg, as a function of the total mass within -R500 (excluding the two galaxies) and 
grouped to a minimum of 8 points per bin. 



double the size of that of iFinoguenov et all is shown in 
Fig.|S] For direct comparison we have derived entropies at an 
enclosed mass of 2.14 X 10^^ /Ito Mo, allowing for our differ- 
ent choice of Hubble constant. Note that in some of the most 
massive systems, the radius enclosing 2.14 x 10^^ h^gMQ 
lies (just) within the cooling region, and hence our derived 
entropy value is effectively extrapolated, using the model 
for the non-cooling component (see section 3). This affects 
six of the clusters which contribute to the highest temper- 
ature bin in the right panel of Fig. |H] and a handful of 
cooler clusters. Our results look significantly different from 



those of IFinoguenov et alJ i2002t) . The entropy appears to 
be non-monotonic, with a minimum value in systems with 
T ~3-4 keV. From a purely phenomenological perspective, 
this non-monotonic behaviour can be under stood in terms 
of mass profiles which take the NFW form llNavarro et alJ 
11997^ 



(3) 



p(r) 



(r/r,)(H-r/r.)2' 



combined with entropy profiles which rise as S oc r^'^, out- 
side some fiatter central core. In groups, the mass shell under 
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consideration here lies well outside Ts, where enclosed mass 
grows roughly linearly with radius, 

M(< r) ~ CrT^l^, (4) 

so that the radius at which an enclosed mass of 2.14 x 
10^"^ h^QplQ is achieved scales as 

r* oc T^,l'\ (5) 
Since the entropy at this radius scales as 
^(r-,) ocTaoori-', (6) 
we expect 

S{T.)^T^,r\ (7) 

which gives a reasonable match to the slope at T < 3 keV, 
as shown in the right hand panel of Fig. |H1 

In more massive systems (M200 > 2 x 10^* h^Q Mq, 
T >3 keV), the mass shell we are studying moves inside 
the scale radius Tb, so that the enclosed mass grows as r^. 
An analysis similar to that above then results in 

5(r.)(xT»or, (8) 

provided that S{r) oc r^'^ at these small radii. However, 
we saw in Fig. Q that for all but the coolest systems, S{r) 
flattens within the cluster core, in which case the positive 
slope of the S'(r«):T relation will be steeper than 0.175, as 
is observed. 

The above arguments seem able to explain the general 
form of the relation in Fig|Hl using just the shape of the 
NFW profile, coupled with entropy profiles similar to those 
seen in simulations incorporating only gravity and shock 
heating. We therefore conclude that there may be no need to 
invoke modification of the entropy profiles by preheating to 
explain this particular result, despite the initially surprising, 
non-monotonic trend. 



5 DISCUSSION 

A substantial amount of theoretical and computational ef- 
fort has been directed towards the problem of similarity- 
breaking in groups and clusters, especially over the past four 
years. We now compare the various models with the results 
above, and with additional information from the companion 
papers to this (Paper I and Paper II), to see how they fare. 

5.1 Cooling models 

A number of author s jBryanl I2OOOI: iPearce et al.lj2000l: 
iMuanwong et al.ll200lt IWu fc Xuell2002bl: iDave et al.ll2002ll 

have suggested that it may be unnecessary to invoke ad- 
ditional heating to explain the entropy floor, since cooling 
will remove low entropy gas from clusters. The counterin- 
tuive fact that cooling can hav e similar effects to heatin g 
in this regard was flrst noted bv lKnight fc PonmanI il997f) . 
who found, on the basis of ID hydrodynamical simulations 
of cluster growth, that cooling acted to flatten the gas den- 
sity proflles, especially in low mass systems, and to steepen 
the L:T relation, but not sufficiently to match observations. 
More sophisticat e d 3D simu lations of c l uster evolution by 
iMuanwong et al.l (I2001I') and iDave et all JioO^) have shown 



that larger effects can be produc ed, dependin g upon the spa- 
tial resolution of the simulations. iDave et ahl show that their 
simulations have converged in this respect, although their 
cooling function does not incorporate emission from metal 
lines, which becomes very significant at T < 2 keV. 

Cooling achieves its effect of breaking the similarity be- 
tween low and high mass clusters, since, at a given gas den- 
sity, the cooling time is shorter at lower temperatures. Hence 
a larger fraction o f the hot baryons cool in groups , compared 
to richer clusters jBrvanl2000t) .l Dave et alJ (12003) show that 
this is able to reproduce the cluster-group L:T relation quite 
well, except that the predicted steepening in the relation 
falls at a slightly lower temperature than that observed. 
They also compar e the predicted gas e ntropies at O.I-R200 
with the results of lPonman et al.1 (|l99flh at T < 4 keV, and 
find reasonable consistency. However, agreement with the 
results from the present study (Figs. ^ and |^, which is su- 
perior to that of IPonman et alJ H1999) in terms of sample 
size and all owance for non-isothermality, is much less good. 
iDave et al] find that S{0.1R2oo) is raised in low temperature 
systems, but converges with the self-similar trend through 
hot clusters for systems with T > 3 keV. In contrast, our 
results show a relation which is flatter t han the self-si milar 
line across the ful l temperature range. iBrvanI (|2Qq3) ^^'^ 
IWu fc Xuel(l2002bh produce results in better agreement with 
our observations, but their analytical model is based on an 
assumed variation in star formation efficiency with system 
mass which is much stronger than that which we derive from 
the subsample of our systems for which we have optical lu- 
minosities (Paper II). 

However, the most fundamental problem faced by 
cooling-only models, as many of their proponents have ac- 
knowledged, is that cosmological simulations without some 
form of effective 'feedback' of energy into the baryonic com- 
ponent as stars form, have a serious 'overcooling' problem, 
which has been recognised for many years. At high redshift, 
the Universe is dense, and a large fraction of the baryons cool 
and form stars wi thin small collapsed haloes. For example, 
IPave et alJ l)2002h find that in the largest systems in their 
simulation (with T ~ 4 keV) almost half of the baryons have 
cooled out of the hot phase, whilst in small groups over 80% 
of the gas has cooled. Unless this cooled gas does not form 
stars, the very high star formation efficiencies implied are 
far above those observed in clusters (c.f. discussion in Paper 
II) or inferr ed globally from the K-band luminosity function 
of galaxies jBalogh et al.ll200lh . Furthermore, recent X-ray 
spectral observations from XMM-Newton and Chandra have 
made it clear that our understanding of cooling gas in the 
local Universe, on scales ranging from clusters to elliptical 
galaxies, is deficient (see e.g. Fabian et al. 200L) . It seems 
likely that cooling rates have been seriously overestimated, 
for reasons which are still under very active debate, but may 
be related to the effects of feedback from AGN, which prob- 
ably exist at the centres of all cooling flows. 

5.2 Preheating models 

Since less energy is required to raise the entropy of gas 
by a given amount when its density is lower, it is ener- 
getically favourable to 'preheat' gas before it is concen- 
trated into a cluster or group potential. A comparison of 
the excess entropy required, with that potentially avail- 
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Figure 8. Left: entropy at the radius enclosing a mass of 3 X 10" Ii^oMq as a function of system temperature (excluding the two 
galaxies). Right: the same data grouped into bins containing 8 points or more, with simple limiting trends, discussed in the text, marked 
for comparison. 



able from supernova explosions associated with star forma- 
tion Ce.g. lValageas fc Silk 199?) shows that a high efficiency 
of heating of the intcrgalactic gas is required even in the 
most favourable circumstances - i.e. only a modest frac- 
tion of the supernova energy must be radiated away. How- 
ever, such high efficiencies have been inferred on the basis 
of observations and mo delling in the case of starburst winds 
JStrickland fc Steven jEooa') . and so are not necessarily un- 
reasonable during the epoch of galaxy formation. 

A large reservoir of energy at high redshift is poten- 
tially ava ilable from the formation of massive bla ck holes in 
galaxies JValaeeas fc Silklll999l: IWu et alJl2000t) . although 
it is at present unclear how much of this energy can be 
coupled into heating of the IGM. Radiative heating from 
AGN cannot achieve the required high entropies, but quasar 
outflows may do so, although the physics and demograph- 
ics involved are still subject to considerable uncertainties 
jNath fc Rovchowdhurll200d) . 

Finally, cosmological simulations suggest that a large 
fraction of the baryon content of the Universe is currently in 
the form of what has been dubbed the 'warm-hot intergalac- 
tic medium', with an entropy of the or d er of 300 keV cm^ 
JCen fc OstrikeJll999l : rDavi et alJliool . lDave et al.l ll200lD 
show that this high entropy gas is primarily heated, not by 
star-formation, but by the effects of shocks generated dur- 
ing the collapse of gas into filaments. It seems unlikely, how- 
ever, that this mechanism can account for the excess entropy 
seen in the inner regions of groups and clusters, since the en- 
tropy of the IGM generated in this way declines steeply with 
redshift JValageas et alJl2002l) . such that it should be well 
below the observed floor value, at the epoch when the gas 
in cluster cores w as accreted. This is c onfirmed by the re- 
cent simulations of lBorgani et al.| ll2002h . for example, which 
find baryon distributions in clusters which are essentially 
self-similar, in the absence of non-gravitational heating pro- 
cesses. 

Pioneering efforts to explore the effects of pre- 
heating through simula ti ons o r analytical treat- 
ments llMetzler fc EvrardI Il994l : iNavarro etall Il995l: 



ICavaliere et al.l Il997^ have been largely superseded by 
more recent work, and we concentrate here on the latter. 
Considerable insight can be gained by analytical models in 
which preheated gas is accreted into clusters. The model of 
Babul, Balogh and coworkers has reached its most advanced 
stage of development in Babul ot al. (2002i). This model 
assumes Bondi accretion llBondiill952f) of preheated gas into 
a potential well represented by a (fixed) NFW profile. This 
accretion is taken to be isentropic in low mass systems, 
with introduction of a shock heated regime, motivated 
by the entropy profiles seen in numerical simulations, for 
systems above some critical mass. The level of preheating 
is tuned to achieve a good match to the L:T relation across 
a wide temperature range (0.3-15 keV). Unfortunately, the 
entropy level required to achieve this (330 keV cm'^), is well 
above that actually observed in the inner regions of groups. 
Moreover, the model predicts large isentropic cores in cool 
systems. The IGM in groups with M200 < 8 x lO^M© 
(corresponding to ~ 2 keV, from our M:T relation in Paper 
I) is entirely isentropic in these models, in strong confiict 
with our results (Fig.0. 

iTozzi fc NormanI J200ll) and lDos Santos fc Dorl i2002l) 

have attempted to evaluate the effects of shock heating on 
the preheated accreting gas. lDos Santos fc Dor j concentrate 
on the scaling properties of the gas entropy immediately in- 
side the accretion shock. They find that this is expected to 
scale with mean system temperature according to (l+T/To), 
where To is an adjustable parameter related to the initial 
adiabat on which the preheated gas lies. By choosing the 
preheating entropy to be 120 keV cm^, corresponding to 
To = 2 keV, they obtain a good m atch to the S{0.1R2oo)'T 
plot of iLlovd-Davies et alj i2000l) . though their curve is 
rather too concave in comparison with the trend seen in 
Fig- El In or der to c alculate the expected L:T relation from 
their model, iDos Santos fc Pore make two key additional 
assumptions: the IGM is assumed to be isothermal, and en- 
tropy profiles, S{r/R2oo), are assumed to have the same 
shape in all systems. Their resulting L:T model steepens 
rather too gently to match the group data (which show a 
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very s teep slope - e.g. L oc T"*'^ from iHelsdon fc PonmarJ 
(|2000|)) entirely satisfactorily. The isothermal assumption 
is in conflict with our data (Fig. |2J and with the results 
of cosmological simulations, which generally show a decline 
in temperature by a bout a factor 3 fro m the inner regions 
to the virial radius (iFrenk et alJll999ll . A picture involv- 
ing self-similar entropy profiles, with a scale factor given by 
(1-1- T/Tq), is in remarkably good agreement with our re- 
sul ts, as shown in Figpi Ifow ever, it should be noted that 
the lDos Santos fc Dord ^|2002^ model predicts only the way 
in which entropy should scale just inside the shock - the 
self-similarity in the radial distribution of entropy inside 
the shock is an assumption, rather than an output of their 
model. 

A more comprehe nsive treatment is provided by 
iTozzi fc NormanI i200ll) . who calculate entropy profiles 
based on solving the shock jump conditions over typical ac- 
cretion histories, for systems spanning a range of final mass. 
They find that shock heating generates entropy profiles with 
a characteristic logarithmic slope S oc r^'^ , outside a central 
isentropic core, which, for a given preheating level, occupies 
a larger fraction of ii2oo in lower mass systems. The slope of 
1.1 agrees well wi th our results (Fig.fH, and with numerical 
simulations (e.g. iBorgani et al.lfeOOlji') . However, as can be 
seen from Fig. we do not see the larger isentropic cores 
in S{r/R2oo) which appear to be a feature of all preheating 
mo dels. 

iNath fc Rovchowdhurvl ll2002h explore the energetics 
and evolutionary history of preheating the IGM through 
outflows from radio loud, and broad absorption line quasars. 
Their conclusion is that AGN could provide the energy in- 
put required to account for the entropy floor. However, in 
the context of our results, their most interesting result is 
that the specific energy injected by AGN into the gas, is 
expected to be higher in low mass systems, for two reasons: 
the incidence of powerful AGN is expected to be higher in 
smaller clusters, and the fraction of their power expended in 
PdV work is larger in smaller systems. This appears to be 
contrary to our entropy scaling results, shown in Figs.|S]and 
|S| which show that at small radii the excess entropy is sim- 
ilar (~ 100 keV cm'^) across a wide range of system masses, 
whilst at larger radii it is actually highest in moderately rich 
(T ~ 3-4 keV) clusters. 

5.3 Star formation models 

A natural development of cooling models, which can help 
to address the overcooling problem, is the incorporation of 
feedback due to star formation. In numerical studies, the 
successful implementation of such a scheme is extremely 
challenging, and has been the Holy Grail of those en- 
gaged in simulations of galaxy formation for some years. 
IVoit fc Bryan (2001) introduced an interesting perspective 
on this issue, in the context of similarity breaking in clus- 
ters, by noting that it makes rather little difference quite 
how effective feedback is in heating the IGM, since any gas 
within virialized systems which has low entropy will have a 
short cooling time, and so is removed from its location near 
the centre of a dark halo either by dropping out of the hot 
phase, or being heated by star formation in its vicinity, such 
that it escapes to large radii. In particular, they noted that 
the entropy of gas with a cooling time equal to the age of the 



Universe is 100 keV cm^, in striking agreement with the 
entropy floor reported in groups and clusters. Since the cool- 
ing time scales as T^^^ /n in systems with T > 2 keV, where 
bremsstrahlung dominates cooling, whilst S = T/n'^''^, it 
follows that a given cooling time is achieved in gas with an 
entropy which scales as 

Scooi oc (^cool^)2''^ (9) 

so that the entropy floor genera ted in this way is ex pected 
to be higher in hotter systems jVoit fc BrvanllioO J) . This 
feature may help to explain the flat trend in entropy with 
temperature apparent in Fig.|S]- note the similarity between 
the dependence of Scooi on T, and the slope seen in Figs. 14151 
and|5| 

T his approach was developed further by IVoit et all 
("2002*), who constructed models in which the entropy distri- 
bution is either truncated below some critical entropy (cor- 
responding to gas cooling out, or being ejected as a result of 
vigorous heating), or shifted by the addition of an entropy 
boost throughout (corresponding to preheating of the entire 
IGM). The two different prescriptions for entropy modiflca- 
tion produce only subtle differences in observable properties: 
shifted entropy models have rather flatter gas density pro- 
files at large radii in poor groups, and also slightly higher gas 
temperatures. These differences are not distinguishable with 
data of the quality used in the present study, and may well be 
too challenging even for future studies with XMM-Newton 
and Chandra. One impressive feature of these models is that 
they contain essentially no adjustable parameters, since the 
entropy threshold is set by equating the cooling time of gas 
to a Hubble time of 15 Gyr. However, the predicted L:T 
relation fails to steepen sufficiently at low temperatures to 
pass through the bulk of the galaxy group points, and the 
M:T relation, whilst steeper than the self-similar relation 
(M oc T^'^), is less steep than that derived from the present 
sample in Paper I. 

One of the most recent numerical studies of cluster 
structure to incorp orate the effect s of both coo ling and feed- 
back is the work of Borg ani et al.l ll200ll. 1200^ . These simu- 
lations explore the effects of setting an entropy floor through 
instantaneous preheating at high redshift, and also through 
heating scaled to the expected supernova energy input in 
overdense regions throughout the evolutionary history of the 
Universe. The effects of radiative cooling are also included 
in one of the supernova heating runs. The main conclusions 
from this work, are that > 1 keV per particle of energy input 
is required to give a reasonable match to the observed steep- 
ening of the L:T relation. The authors point out that such a 
large energy boost appears to exceed what is expected from 
supernova input alone, even if a high heating effeciency of 
the IGM is assume d. 

iBorgani et al.l flnd that a somewhat larger suppression 
of Lx in groups is achieved by injecting the same amount 
of energy progressively, compared to global preheating of 
the IGM. Entropy pro files generally agree well with the 
iTozzi fc NormanI (1200 ll) slope of 1.1, but are flattened in 
low mass systems by the effects of energy injection, as are 
gas density proflles. The gas temperature proflle becomes 
more centrally peaked in low mass haloes subject to ad- 
ditional heating. The M:T relation is found to be little af- 
fected by heating, either in slope or normalization. However, 
cooling is found to decrease the normalization and steepen 
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the slope of the relation, bringing it into better agreement 
with observations. In comparison to these results, we cer- 
tainly observe flatter gas density profiles in cool systems 
(Paper I), and Fig. |5| provides some tentative evidence that 
the entropy slope may be rather shallower in lower temper- 
ature (T < 3 keV) systern s. In the case of t he M :T relation, 
iLlovd-Davies et alJ i2002l) and lVoit et"al] J2002ll show that 
systematic variations in the concentration parameter of dark 
matter haloes with mass, may play a key role in steepening 
its slope, as discussed in Paper I. 



6 CONCLUSIONS AND SUGGESTIONS 
6.1 Existing models 

Drawing the above results together, it seems that cooling- 
only models cannot provide a viable explanation for the 
similarity-breaking observed in clusters, unless one is will- 
ing to admit the presence of very large quantities of baryonic 
dark matter, which would have to dominate the baryon bud- 
get in groups. 

Preheating models appear to have the generic property 
of generating large isentropic cores in low mass systems, 
in conflict with our observations. This seems to be an in- 
escapable feature of simple preheating models in which the 
IGM is raised uniformly and 'instantaneously' to a high adi- 
abat, since such gas can only be shocked when falling into 
potential wells deep enough for its motion to be supersonic. 
More complex preheating models may enable this problem 
to be circumvented, as we discuss below. 

Models involving a mixture of cooling and star forma- 
tion (or possibly AGN heating) appear more natural and 
promising. However, a number of features of our data ap- 
pear to conflict with all models proposed to date: 

(i) The very large entropy excesses seen at large radii 
(Figs. |S| and [7J are a suprise. Models tend to show entropy 
enhancement in the inner regions of low mass systems, with 
a normal shock-generated entropy profile re-establishing it- 
self at larger radii. 

(ii) Closely related to this, is the fact that the entropy 
profiles appear to be approximately self-similar apart from 
a normalization constant, and in particular, that larger isen- 
tropic cores are not seen in galaxy groups. In fact. Fig. 
shows that the lowest temperature systems actually have 
entropy profiles which appear to drop all the way into the 
centre, u nlike hotter systems. The fa ct that the scaling sug- 
gested bv lDos Santos fc Dorel ll2002l) brings our profiles into 
good agreement (Fig. |2J is a puzzle, given that this scaling 
is really not justified by their model as a scale factor for 
entropy profiles. 

(iii) The temperature profiles shown in Fig. |21 are also 
not quite what is expected from the models. Any mechanism 
which gives an entropy boost should produce the strongest 
results in low mass systems. Since a rise in entropy has to 
be coupled with the maintenance of hydrostatic equilibrium, 
the natural consequence is a rise in central temperature, cou- 
pled with a decrease in density. The density drop is certainly 
observed, and there are indications that the temperature has 
been raised in cool systems outside the core, but the tem- 
perature profiles in groups seem to lack the expected central 



cusp. This is related to the lack of an entropy core in groups 
referred to in point (ii), above. 

(iv) The behaviour of the gas core radius, discussed in 
Paper II, is similarly unexpected. As a fraction of -R200, the 
gas core radius is approximately constant at ~10% over most 
of our temperature range, but falls dramatically, by an order 
of magnitude, in group s. The only study to predict this sort 
of behaviour is that of IWu fc Xul ll2002ari . who show that 
fitting a convex profile which steepens progressively, with a 
beta model, can lead to smaller fitted core radii as the outer 
radius of detection shrinks. In Paper I we reported tests 
with fits to real data for two clusters, truncated at different 
radii, which suggested that this effect is not dominating our 
fits. However, X-ray surface brightness profiles extending to 
larger fractions of -R200 for groups are required to definitively 
settle this issue. 

6.2 New possibilities 

What might these disagreemen t s with the models be point- 
ing towards? iFinoeuenov et all (l2002l) argued that the large 
excess entropy at -R500 in groups indicates that the entropy 
profiles in cool systems are dominated throughout by the ef- 
fects of preheating. Two features of our results make us un- 
easy about this conclusion. Firstly, these excesses are now 
seen (Fig. |7J to extend up to moderately rich clusters, and 
to involve entropies se veral times those expec ted in the IGM 
at the present epoch JValaeeas fc Silk^ll999^ . Secondly, the 
slope of all our entropy profiles (Fig. seem to scatter 
about the value predicted from shock heating. This is an 
uncomfortable coincidence if shock heating plays no role in 
establishing them. It is certainly possible to devise preheat- 
ing scenarios which could reproduce any of our profiles, by 
varying the history of energy injection, and the density gra- 
dient of the gas into which this energy was injected, in the 
vicinity of systems of a given mass. However, such a solution 
would require a large amount of fine tuning. 

If, on the contrary, we wish to retain shock heating as 
the basic mechanism responsible for generating the rising 
entropy profiles seen outside the core in both clusters and 
groups, how can this be achieved? The post-shock temper- 
ature is essentially determined by the infall velocity of gas 
into a halo, which could be reduced by pressure forces, but 
not increased above the free-fall value. Hence, the only way 
to raise the whole entropy profile in low mass systems, in 
the way implicit in Fig. Q seems to be to reduce the den- 
sity of the pre-shock gas accreting into lower mass systems, 
relative to that falling into rich clusters. In simple spheri- 
cal collapse models, the matter turning around and accreting 
into haloes at a given epoch is expected to have a given over- 
density. However, in reality, cosmological simulations show 
us that accretion is far from spherical, and that much of it 
takes place along filaments. 

This change in topology does not in itself change the 
relationship between the mass of virialized systems and the 
density of the gas accreting into them, since structure for- 
mation is (almost) self-similar. Bigger systems have bigger 
filaments, and the mean density of the accretion flow at a 
given epoch is still set by the overdensity at turnaround, 
which is independent of system mass. However, this situa- 
tion changes if preheating of gas in the filaments introduces 
an additional scale into the problem. In general, the scale 
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height of gas in filaments will be set by hydrostatic equi- 
librium of the ga s in the local gravitat ional potential of the 
dark matter fe.g. lValageas et al.ll2002ll . The temperature to 
which gas is shock heated during the collapse of filaments 
will scale with their mass, and hence the effect of any ther- 
mal input into the IGM from galaxy formation or the growth 
of supermassive black holes will be more prominent in the 
smaller filaments associated with lower mass systems. Pre- 
heating will act to increase the scale height, and hence to 
reduce the density of gas in filaments. Even in filamentary 
structures which are still collapsing, preheating of the gas 
will increase its pressure, and retard its collapse relative to 
that of the dark matter. 

The calculations of iTozzi fc NormarJ (1200 ll) show (al- 
beit in the spherically symmetric case) that even in the pre- 
heated case, the accretion shock rapidly becomes strong once 
it is established, so that the post-shock entropy of gas ac- 
creting into a cluster will scale as T200 /n^^^ , where ni is the 
pre-shock density of the gas, and T200 is the virial tempera- 
ture of the cluster. Hence, if preheating reduces n\ in lower 
mass systems, relative to rich clusters, then the entropy will 
scale sub-linearly with T200, which is what we observe. 

It is worth noting some features of this model. Firstly, 
the effect we invoke is intimately linked to the fact that 
the accreting gas is largely confined to filaments, and hence 
will not be seen in analytical a nd semi-analytical models, 
such a s those of iTozzi fc NorniTan (^2001) and Ba bul et alJ 
i2002h . which assume spherically symmetric accretion. In 
the spherical case, the gas has nowhere to expand to in re- 
sponse to the initial preheating. Consequently any entropy 
rise due to preheating serves to raise the temperature of the 
gas, rather than to lower its density. In the lowest mass sys- 
tems, this temperature rise can delay shock formation, and 
hence result in an isentropic core, but in high mass systems 
(and at large radii in low mass ones) it has essentially no ef- 
fect, since the pos t-shock temperature is just T = ^fimpvf 
l|Tozzi fc Normanll200ih . where Vi is the velocity at which 
gas, with mean particle mass fimp, fiows into the shock. In 
contrast, if preheating serves primarily to reduce the den- 
sity of the pre-shock gas, then this will lead to a rise in 
post-shock entropy whatever the shock strength, since (for 
a given shock Mach number) the post-shock density simply 

— 2/3 

scales with ni , so that the post-shock entropy S2 (x n-^ 

Secondly, this mechanism is a very efficient way of rais- 
ing the entropy of int racluster gas for two reasons. As has 
been noted previously jPonman et al .ll999^ a given injection 
of energy produces a larger rise in entropy when the gas den- 
sity is lower. Hence it is more effective to deposit energy into 
low overdensity structures such as filaments, than into fully 
collapsed systems like clusters. Moreover, if this injection 
places the gas onto a higher adiabat by reducing its den- 
sity, then the effects of the accretion shock serve to raise it 
further. For example, Dos Santos & Dorc (2002) show that 
the ratio between the post- and pre-shock adiabats is well 
approximated by 

+ ™ 

where A is a constant of order unity. So, a rise in S'l by some 
factor, due to a density change (and hence with no change 
in Ti) driven by preheating, will boost S2 by the same fac- 
tor. Hence the shock has a multiplier effect, and if S2/S1 is 



large, a modest rise (in absolute terms) in Si, may result 
in a much larger rise in 82- Since the entropy will typically 
jump by an order of magnitude in the accretion shock, as its 
temperature is raised from ~ 10^ K to 10^-10* K, entropy 
excesses of the magnitude reported here within clusters (i.e. 
100-1000 keV cm^), might be generated from a rise in en- 
tropy of only '--^ 10- 100 keV cm'^ within the filaments which 
feed them. 

Thirdly, if preheating acts primarily to reduce the den- 
sity of the pre-shock gas, rather than to increase its tem- 
perature, then the Mach number of accretion shocks is not 
much reduced, compared to the unhealed case. The post- 
shock entropy profiles are therefore almost entirely due to 
the evolution of the accretion shock, and S{r) cc r^'^ pro- 
files are expected, outside a small unshocked central core, in 
good agreement with observations. 

Whilst the viability of this model clearly cannot be in- 
vestigated by spherically symmetric analytical treatments, 
one might hope to find evidence to test it from 3-dimensional 
simulations of cluster formation. Rather few studies have 
been published which incorporate the effects of feedback 
into the IGM, together with a study of its effects on the 
entropy profiles of cluster s. T he results from the wo rk of 
iBorgani et all (I200lll2002l') and lMuanwong et alj ll2002h . are 
not very encouraging. In both these studies, the effects of 
heating of the IGM serve to fiatten the entropy profiles in 
clusters, with the effects being more pronounced in lower 
mass systems. This accords better with observations than 
do most 1-dimensional analytical models, insofar as distinct 
isentropic cores are not generally present. However, the pro- 
gressive fiattening of S{r) towards lower mass systems does 
not agree with our observed properties (Fig. 0. 

It is premature to conclude at this stage that the 
preheated filament model can be rejected. Feedback is 
very poorly understood, and its incorporation into nu- 
merical codes is notoriously difficult. In the study of 
iMuanwong et all (l2002h . feedback was included by sim- 
ply raising the tem perature of all gas pa rticles by 1.5 keV 
at a redshift of 4. iBorgani et al.l (1200 adopted the ap- 
proach of imposing an entropy floor on the gas in overdense 
(5 > 5) regions at a variety of redshifts from 1 to 5, whilst 
IBorgani et alJ ll2002h attempted to model the effects of en- 
ergy injection from galaxies in a more realistic way by using 
a semi-analytic scheme to predict the supernova rate, and 
sharing the resulting energy amongst gas particles in regions 
with S > 50. Unfortunately, the interplay of feedback and 
cooling is extremely complex, and hard to model. Results 
tend to depend strongly on spatial resolution, and current 
simulations (including those of Borgani et al. and Muan- 
wong et al. ) generally fail to reproduce the observed low 
cooled fraction of baryons. A wider exploration of the effects 
of different feedback schemes, and a study of the density 
of gas flowing into accretion shocks in simulations, should 
clarify whether the mechanism proposed here is actually the 
dominant effect in modifying the entropy of intracluster gas. 
At the same time, higher quality observations from XMM- 
Newton and Chandra will provide more robust determina- 
tions of observed entropy profiles in a wide range of virial- 
ized systems. Initi al results from XMM-Newton for a small 
nurnber of groups dMushotzkv et al]l2003l : iPratt fc ArnaudI 
I^OSD confirm that no central isentropic core is present in 
these systems. 
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